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Abstract 

In this paper, the effect of the control on bipartite entanglement is discussed 
from a geometric viewpoint for a nuclear magnetic resonance (NMR) system 
as a model of the n-qubit control system. The Hamiltonian of the model is 
the sum of the drift and control Hamiltonians, each of which describes the in- 
teraction between pairs of qubits (or ^-spins) and between one of qubits and 
an external magnetic field, respectively. According to the bipartite partition 
(C2)®" = (C^)®^ (g) (C^)®"^ with ^ + m = n, the Schrodinger equation for the 
NMR system is put in the matrix form. This paper gives a solution to the 
Schrodinger equation with the assumptions of small coupling among qubits and 
of constant controls. The solution is put in the form of power series in small 
parameters. In particular, in the case of the two-qubit NMR system, the drift 
and control Hamiltonians are shown to be coupled to work for entanglement 
promotion, by examining solutions to the Schrodinger equation in detail. The 
concurrence, a measure of entanglement, is evaluated along the solution for a 
small time interval in order to observe that the control effect appears, not at 
the first-order terms in t, but at the higher-order terms in t. The evaluated 
concurrence also suggests which control makes the two-qubit more entangled or 
less. 
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1 Introduction 

Entanglement plays many roles in quantum computation and quantum information 
theory, and is a resource for quantum communication. Recently, the geometry of 
entanglement has been developed [9,11-13,18-22]. The algebraic study of entanglement 
also has been made to seek for invariants that quantify multiqubit entanglement [4,10, 
26]. The geometry of entanglement in few-qubit systems was studied in [18,22] and 
in [20,21] in terms of the Hopf fibration and the complex projective spaces, respectively. 
The first author [11] showed that the concurrence for the two-qubit system can be 
characterized as the coordinate of the factor space of the state space by the action 
of the local transformation group. He continued to study bipartite entanglement in 
multiqubit systems [12]. In particular, the sets of separable and maximally entangled 
states in bipartite entanglement are identified, and the concurrence is extended to 
a measure of bipartite entanglement in multiqubit systems. Further, the distance 
between the set of states of prescribed entanglement and the set of separable states 
is measured. The bipartite entanglement is measured also for the process of Grover's 
search algorithm [13]. For the bipartite entanglement, see also [9,19]. 

On the other hand, from a practical standpoint, it is of much importance to control 
quantum states in order to create entangled states. Many of papers on control problems 
for multiqubit systems deal with control systems on unitary groups [1-3,5,6,15-17,23]. 
The theory of control on compact Lie groups has been estabhshed [14] . A sufficient con- 
dition for controUabihty of the control system on SU(A'") was given in [3]. The Cartan 
decomposition of the Lie algebra and the symmetric space are techniques frequently 
used for studying control systems [15-17]. 

Two notions as to controUabihty of multiqubit systems were defined and studied by 
Albertini and D'Alessandro [1,2], which are "pure state controllability" and "equivalent 
state controllability". In [2], equivalent state controllability was shown to be equivalent 
to pure state controllability. In the case of multiqubit systems lying in electro-magnetic 
fields, they [1] gave a necessary and sufficient condition for the system to be (pure state) 
controllable. The condition is described in terms of the graph associated with the 
interactions of qubits, which is related with the Lie algebraic structure of the system. 
Their method is based on the control theory on compact Lie groups [6, 14]. As will 
be pointed out in Appendix B, their model Hamiltonian is a bit different from that 
adopted in this article. 
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The aim of this article is to study entanglement and control of nuclear magnetic 
resonance (NMR) multiqubit systems from a geometric viewpoint. Since the notion of 
entanglement cannot be defined in terms of the whole transformation group, the control 
theory on compact Lie groups will not fit in with the study of entanglement and control. 
The state space of a finite-dimensional quantum system is a finite complex vector space 
with a constraint required for the probability density. In this article, the state space 
is a linear space of complex matrices with the constraint from the probability density. 
The reason for the choice of the present state space is that the entanglement measure 
like the concurrence is easy to treat on this state space. The state space admits the 
two-sided action of the local transformation group U(2^) (g) U(2™), where m and i are 
determined by the bipartite partition of the multiqubit (C^)®" = (C^)®'^ O (C^)®"^. 
As will be stated later, the entanglement of the multiqubit system is invariant by the 
action of the local transformation group. Hence, in order to make the multiqubit more 
entangled, the control is required to move the state in the direction transverse to the 
orbit of the local transformation group. However, for NMR systems to be treated in 
this article, the control vectors generated by the control Hamiltonian are tangent to 
the orbit mentioned above, so that they would not make the system more entangled, 
if no coupling were made with other vector fields. Since the total Hamiltonian of the 
NMR system is the sum of the control and drift Hamiltonians, one may expect the 
coupling between them to occur. The objectives in this article are to set up the control 
problem on the above-mentioned state space and to show, by solving the Shrodinger 
equation with constant controls, that the control and drift Hamiltonian are coupled 
together to give rise to vector fields transverse to the group orbit for the entanglement 
promotion of the multiqubit system. 

The organization of this paper is as follows: Section 2 contains a geometric setting 
up of n-qubit systems. According to the bipartite partition (C^)®" = (C^)®^ <S> (C^)®"^ 
with i+m, — n, the state space M is viewed as the space of 2^ x 2"* normahzed matrices 
C, on which the local transformation group G — U(2^) ® U(2"*) acts. The action of 
G leaves invariant the bipartite entanglement of the n-qubit. Linear operators on the 
n-qubit system are mapped to vector fields on the state space M. Lie brackets among 
those vector fields are given explicitly. In Section 3, the Hamiltonian of an NMR 
model is given, which is split into two operators, the drift Hamiltonian Hd and control 
Hamiltonian He- The Hd and He describe the interactions between pairs of n qubits 
and between one of qubits and an external magnetic field, respectively. In addition. 
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a solution to n-qubit NMR systems with a certain type of drift Hamiltonians is given 
with the assumption that the couphng constants between qubits are small enough. The 
solution is put in the form of power series in small parameters related to the coupling 
constants. It will be shown that the vector fields associated with the drift and control 
Hamiltonians are coupled actually. In Section 4, the coupling between those vector 
fields is extensively studied on the two-qubit system to observe that entanglement 
is actually promoted. Discussion runs as follows: A solution C{t) of the Schrodinger 
equation in the matrix form is obtained as a power series in the small coupling constant. 
The solution is examined in detail to show that the vector fields associated with the 
control and drift Hamiltonians are coupled indeed to induce vector fields transverse to 
the G-orbit. The effect of the control on entanglement is verified by estimating the 
difference between the concurrences of the solutions for the controlled and uncontrolled 
systems for a small time interval in the case that the initial state is a diagonal matrix 
A = diag (Ai, A2). The result is that the effect of the control is so slow that it emerges 
at the third-order term in t. Concluding remarks are made in Section 5, in which the 
Cartan decomposition method is related to the present method. Appendix A gives a 
list of Lie brackets among vector fields related with the drift and control Hamiltonians 
for the two-qubit NMR system. Appendix B contains the proof of the controllability 
of the n-qubit NMR system. The method for proof is different from that in [1] . 

2 Geometric setting for bipartite entanglement 

2.1 The state space 

We start with geometric setting for the n-qubit system after [12]. The Hilbert space 
for the n-qubit system is (C^)®" = (g) • — (g) C^^ , whose state vectors 1^^) are put in 

n times 

the form 

J2 Cj\ji)^---^\jn) with ^|cj|2 = l, 

J=Oi,- Jn)e{o,i}" J 
where {|0), |1)} denotes the computational basis of a single qubit system C^. We 
denote by {| J)}jg{o,i}" the standard orthonormal basis |ji) |j„) of (C^)^". The 

Hermitian inner product of l^') = ^jCj|J) and |$) = '^jdj\J) is then given by 

According to a bipartite partition of the n qubits, the Hilbert space (C^)'^'^ is 
decomposed into the tensor product (C^)®^ (C^)®'" with i + m, = n. We here assume 
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that < i < m without loss of generahty. A state vector |^) of the n-qubit system 

is separable with respect to the bipartite decomposition (C^)®" = (C^)*^^ (C^)®™, if 
there are two unit vectors E (C^)®^ and l^a) e (C^)®"* such that |^') = |^'i)(g)|^'2)- 
A vector is said to be entangled if it is not separable. 

According to the bipartite decomposition (C^)®"" = (C^)®^ (C^)®"^, a state vector 
1^) can be rewritten in the form 

where A = (ji, • • • ,je) and B = (j^+i, ■ • ■ This expression gives rise to the 

isomorphism of (C^)®" to the vector space C^*^^" of 2^ x 2"^ complex matrices, 

, : (c^)®- C''-^- ; 1^) = J]c^s|A) ® |S) ^ C = (cab). (2.1) 

A,B 

The C^^^^"' is endowed with the inner product given by tr(C*C2) for Ci, C2 G C^*^^"". 

Since |^) is normahzed, the corresponding matrix C = '-(l^)) is subject to the 
condition tr(C*C) = 1. Thus, the state space for the n-qubit system is defined to be 

M := {C e C^'^^^"" I tr(C*C) = 1}, 

which is diffeomorphic with the (2""'"^ — l)-dimensional sphere. The state space M 
becomes a Riemannian manifold equipped with the standard metric 

{X,Y)c := ^tr{X*Y + Y*X), X,Y e TcM, (2.2) 

where the tangent space to M at C is identified with 

TcM ^{X e C^'^^" I tT{X*C + C*X) = 0}. 

2.2 Bipartite entanglement 

In this subsection, we make a review of a measure of bipartite entanglement introduced 
by the first author [11,12]. 

The separability condition of a state vector |^) can be expressed as rank CC* — 1, 
where C is the matrix corresponding to |^). Since the 2^ x 2^ matrix CC* is positive 
semi-definite together with tr (CC*) = 1 and I <m,, the condition of rank CC* = 1 is 
equivalent to det(/2< — CC*) = 0, where denotes the 2^ x 2^ identity matrix. The 
first author [11, 12] showed that the quantity 

F(C) := det(/2^ - CC*) (2.3) 
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serves as a measure of bipartite entanglement, and as an extension of the concurrence 
which is well-known as a measure of entanglement in two-qubit system. In fact, one has 
F{C) = det(CC*), the square of the concurrence, for the two-qubit with i = m = 1. 
The quantity F{C) is non-negative, and vanishes if and only if C is separable. Further, 
the F takes the maximal value (1 — 1/2^)^* if and only if all eigenvalues of CC* coincide. 
If F attains the maximal value at C, the matrix C G M or the corresponding state 
vector |\&) is said to be maximally entangled with respect to the bipartite decomposition 

2.3 U(2^) U(2'^)-action 

The action of the group G = U(2^) ® U(2'^) on M C C^'^^™ jg defined by 

C ^ gCh^, (2.4) 

where g G U(2^) and h G 11(2™). This action is isometric with respect to (12.21) . We 
remark here that the map l is G-equivariant, namely, L^^{gCh^) = {g ® h)L^^{C). As 
is easily verified, the function F is invariant under the G-action, so that the G-action 
(12.41) does not change bipartite entanglement. 

We now consider the quotient space M/G. An arbitrary matrix G G M is decom- 
posed into the product 

G = g{k, 0)h^, g G U(20, h G U(2™), A = diag (Ai, ■ ■ ■ , A^.) , (2.5) 

where is the 2^ x (2"* — 2^) zero matrix, and Af , ■ ■ ■ ,X'^e are the eigenvalues of GG*. 
In the following, we assume that Ai > • ■ ■ > A2« > 0. The quantities Ai,--- , A2^ 
are called the singular values of G. The expression (12. 5p is equivalent to the Schmidt 
decomposition of a state vector |^). From this decomposition, the quotient space M/G 
can be identified with 

|(Ai, • • ■ , A2O G M^' I Ai > ■ ■ ■ > Aa^ > 0, A? + • ■ ■ + A^, = 1 1 . 

The sets of separable and maximally entangled states project through the natural pro- 
jection M ^ M/G ; Gh^ (Ai,-- - ,A20tothepoints(l,0, -- - , 0) and (1/v^, - - - ,1/^2^) G 
respectively, both of which lie on the boundary of the quotient M/G. 
The state space M is stratified into strata, according to G-orbit types [12], which 
are determined by the isotropy subgroups Gc ■= {g®h G G | gGh^ = C}. The isotropy 
subgroup and the orbits are already studied for four-qubit systems in detail, and for 
multiqubit systems in brief [12]. We now summarize the results in a refined form. 
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Proposition 2.1. Let N be the number of distinct singular values of C, and mj the 

multiplicity of its j-th largest singular value. 

(1) If CC* is non-singular, the isotropy subgroup Gc is isomorphic to 

(U(mi) X ■ ■ • X U(miv)) ® U(2" - 2^), (2.6) 

where the product U(mi) x ■■■ x U(m7v) is viewed as a subgroup o/ U(2^), and the 
G-orbit Oc '■= {{g ® h)G\g ^ h E G} is diffeomorphic to 

G/Gc ^ u{r) ^ {U{mi)x---xU(mjv)) 

VAC'"'), (2.7) 

where V2^(C^'") denotes the Stiefel manifold of orthonormal 2^ frames in C^™. 

(1-1) In the generic case that singular values ofG are all distinct and non-zero, Gc 
and Oc take, respectively, the form 

= T'' ® U(2" - 2^ and Oc ^ lii2^) x^,, V2iiC'^). 

(1-2) For a maximally entangled state C, Gc and Oc become, respectively, 

Gc = U(2^) ® U(2'" - 2^) and Oc ~ ^HC'"). 

(2) IfGG* is singular, the group Gc is isomorphic to 

(U(mi) X ■ ■ ■ X U(m^„i) x U(m^)) ® U(2"^ - 2^ + m^v), (2.8) 
and the orbit Oc is diffeomorphic to 

G/Gc ~ V'2^-m^r(C^ ) X(U(mi)x---xU(mjv-i)) ^2^-miv(C^ )• (2-9) 

(2-1) In particular, if G is separable, Gc and Oc are, respectively, of the form 
Gc = V{2' - 1) X U(2™ - 1), and Oc ~ S^'^"-^ Xu(i) S^^^'-^ 

2.4 Vertical and horizontal subspaces 

This subsection deals with the infinitesimal action of G = U(2^) 011(2'"). For ^ G u(2^) 
and rj G u(2™), the fundamental vector field associated with ^ ® /2m + /2'; ® ''7 G = 
Lie(G) is determined through 



e*«Ge*'^ =^G + Gr/^. (2.10) 

t=o 



^The symbols = and w are used to denote the isomorphism and the difFeomorphism, respectively. 
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The vertical subspace Vc is defined to be the tangent space TcOc at C to the orbit 
Oc, and formed by the fundamental vector fields evaluated at C, 



Vc = {X^^i,^ + r,^,iC) = + Cr^^ I e G u(20, G u(2"^)}. 



The dimension of Vc is given by dim Vc = 2^^ + 2^™ — 1 — dimG^. In particular, from 



Proposition 12.11 we have 

f2^+i + 2"+i-3, if C is separable, 
dimVc = < , „ (2-11) 

I 2*^+™+^ — 2 , if C is maximally entangled. 

The horizontal subspace He of TcM is defined to be the orthogonal complement 
of Vc with respect to the metric (12. 2p . A vector X G TcM is horizontal, if and only if 
{X,C,C + C?7^) = for any ^,ri E u(2^), so that the horizontal subspace at C G M is 
expressed as 

Hc = {X e TcM I XC* - CX* = 0, C*X - X*C = 0}. 
From (12. lip together with dimif,^ = 2^+™+^ — 1 — dim VJ^, we have 

dim He = 



2(2^-l)(2™-l), if C is separable, 
2^^ — 1, if C is maximally entangled. 



2.5 The Lie algebra of linear vector fields on M 

In discussing control problems, we need to calculate Lie brackets of vector fields on 
C^'^^^™. In particular, we have to work with the Lie brackets of vector fields of the 
form 

Xa^b{C):=ACB^, Cee'^'^, 
where A G and B G C2'">^2- 

need to be Hermitian or skew Hermitian at 
present. The Lie bracket of ® 5i, ® ^2 e C^''''^' ® C2">^2- defined to be 

[Ai ® 5i, A2 ® B2] := [Ai, A2] ® B^B2 + ® B2]. (2.12) 

Since this bracket is verified to satisfy the Jacobi identity, the linear space C (S> 
(j-^2'"x2'" endowed with a Lie algebraic structure. This Lie algebra is isomorphic to 
0[(2",C) through 

/ auB ■■■ 

f : C^'x^' (^C^-'x^"' ^C^"^^"- ■ A0B ^ : : 

ya2« iB ■ ■ ■ 2^-8 / 
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In fact, / is shown to be bijective and homomorphic, f{[Ai ® Bi, A2 ® -B2]) = [/(^i ® 
B,)J{A2®B2)]. 

Now we are in a position to state the following proposition: 

Proposition 2.2. For Ai ® Bi,A2 ® B2 e C^'^^^ ® C^'"^^^-^ ^/^^ j^^^ bracket of the 
associated vector fields Xa^^Bi o,nd Xa2®B2 on C^*^^™ is given by 

[Xai(^Bi,Xa2»B2] = —XlAi»Bi,A2<2>B2]- 

This means that the correspondence A^ B \—>- Xa^b gives rise to an anti-isomorphism 
of (C^' ® C^", [•,•]) to the Lie algebra {Xa.^b\A G C^''^'^\B e C2'">^2-| ^^^^^^ 
vector fields on C^'^^^™'. 

Corollary 2.3. /n t/ie two-qubit system with i = m = 1, one has 

for j, j', k, k' = 0,1, 2, 3. Here the aj,j = 1, 2, 3, denote the Pauli matrices defined as 
1 /O l\ 1 / O A 1 A 



■ 2 VI oy ' ■ 2 v-« oy ' • 2 VO -1 



Proof of Proposition \2.2\ . Let v^^^^ denote the flow generated by the vector field 
Xa^b on C^'x^^. We identify (C^)®" with the space C^" of 2"-dimensional column 



vectors, and the isomorphism t : (C' 



c 



with an isomorphism of C^" to 



C 



2'^x2'" 



We note here that the following diagram commutes, 

^2^x2™ action of AiXi-B (j-^2*x2™ 



This means that 



action of f{A<S$B) 



c 



2" 



{l ^)*Xa®b(C) = ^(C)) and t 



,tfiA<g>B) Q 



(2.14) 



where F/(a®b)(c) = f{A (g) i?)c with c G C^". 
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By the definition of tlie Lie bracket togetlier witli fl2.14l) . we obtain 



dtds 



t=s=0 



i O (e"''^(^2«5B2) o g-t/{Ai®Bi) Q ^sf{A2(X)B2) o gt/(Ai®Bi)^ ^ t"^(C) 



t=s=0 



dtds 



Furtlier, by the formula [1^=^,1^=2] = ~^[Si,H2] ^'^Y '='i)'='2 G C , one has 

= '•* [^/{Ai^Bi),'i^/(A2»B2)](''""^(C)) = (-'^[/(Ai®Bi),/{A2®B2)])(''""^(C)) 
= — V}-([A^^B^^yi2^B2])('- ^(C)) = —X[Ai^Bi,A2(^B2]iC). 



2" 



Since the map A ® i? 1— > F/ (A(g)_B) is bijective to the space of hnear vector fields on C 
so is the map B Xa(^b to the space of linear vector fields on C^^^^"". This ends 
the proof. □ 

We now apply 02.131] for the two-qubit system. We take, for instance, Ai ^ Bi = 
ia^ C?) cr^ and A2 ^ B2 = icri ® I. By Corollary 12.31 the Lie bracket of the associated 
linear vector fields is given by [Xi^.^^^.^, Xi^^^g^i] = X[^^^^^]^^^ + X^^^^^^^^j^ = -Xi„^^„^. 
Other Lie brackets will be found in Appendix A. 

3 Controls in NMR systems 

Sections [3] and H] deal with control problems in n- and 2-qubit systems, respectively. 
In this paper, we adopt a nuclear magnetic resonance (NMR) system as a quantum 
computation model from n-qubit systems [16]. 

The total Hamiltonian operator H of the NMR system is the sum of two operators, 

H := Hd + He, 

where Hd and He are called the drift and control Hamiltonians, and expressed as 

H,:= Yl (3.1a) 

l<o</3<n 

TT \ ^ / («) (a,n) . (a) {a,n)\ 

Hc-= 2_^\v\'a\ '+t;^V^ M, (3.1b) 

a=l ^ ^ 
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respectively. Here, Ja/s G IR are coupling constants determining the strength of the 
interactions between the a-th and /3-th qubits. The : R ^ R are (time- 

dependent) controls acting on the a-th qubit. The crj"'"^ are defined as 

(3.2) 




for J = 1,2,3, and a — 1, • • • ,n, with / the 2x2 identity matrix. The controls 
and are assumed to be piecewise constant functions in t. 
The Schrodinger equation for the NMR system is expressed as 

d 



dt 



I*) = -iH\-^), (3.3) 



where the natural unit system has been adopted, so that the Plank constant h is 
equal to one. Since we are interested in bipartite entanglement, we have to put the 
Schrodinger equation in the form of matrix equation. According to the isomorphism 
^(^2j®n ^ £2'x2"' ^-^j^ n — £ + m,, the vector fields associated with the drift and control 
Hamiltonians are expressed, respectively, as 



_.(a,m) (/3,m)ir 
^3 J 



i<a<t3<e e+i<a<p<e+m 
- ^ C [af -^'™)]^, (3.4a) 



e+i</3<e+m 



-iH,{C) - t^®^""'^ ®^<-® ^^^'""^] ^ + H ^ <g, <g, J®(m-/3+£)|r^ 

l<a<i e+i<f3<e+m 

(3.4b) 

where 

and where the symbols with square brackets like [(Tg"'^^] denote the Kroneckcr products 
of matrices corresponding to the tensor products concerned. The Schrodinger equation 
on M is then put in the form 

^^-tH,{C)-zH,{C). (3.5) 

As is easily verified, —iH^ and —iHc are vector fields on M. In what follows, we will 
refer to —iHa{C) and —iHc{C) as the drift and control vector fields, respectively. 
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3.1 A solution to the NMR system 

We show that Eq. (13.51) can be solved with the assumption that the couphng constants 
are small enough and that the controls are constant. In view of (13.41) . we may put the 
drift and control vector fields in the form, 

N 

-iH,{C) = -iY. e,J:f^C{J:i'Y, (3.6a) 

k=l 

-iH,{C)=iiC + Cil, (3.6b) 

respectively, where e/c, /c = 1, ■ ■ ■ , A^, are small parameters, where G C^*^^™, Sg'^^ G 
(^2^x2™^^ _ . . . ^ constant Hermitian matrices, and where .^i G su(2^),i^2 G 

su(2'^) are constant controls. The Schrodinger equation we are to solve is 



k=l 

N 

= - y^.^kX.y(k)^y(_k) (C) +X^j0/2„(C) +X/^^0^2(C). (3.7) 

k=l 

Since the coupling constants Ek are sufficiently small, we may suppose that the solution 
C(t) to (13.71) can be expanded into a power series in Sk, k = 1, ■ ■ ■ , N, 



C{t) = J2 ^'^Cr.it) with := e^l' ■ ■ -e^. 

ri=(ni,---,njv)e{Z>o)^ 



By substituting this expansion into (13.71) and comparing the left and right hand-sides, 
we obtain the series of differential equations 

C^ = -t ^?Cr,^.,{4'Y + ^iCr^ + C^^I for n^O, (3.8) 

k=l,--- ,N, 
nfe^O 

where ei = (1, 0, ■ ■ ■ , 0), ■ ■ ■ , e^v = (0, ■ ■ ■ , 0, 1) G (Z>o)^. 

To solve the above equations, we introduce new unknown functions by Qn{t) = 
(e~*^^ ® e~*^2) ■ Cn(t). Then Qnif) with n G (Z>o)^ prove to satisfy the equations 

Qo = 0, 

fc=l,--- ,N, 

= -i ^ Adg-tei^g-t52(sf^ ® S2^^)(5'n,-efe for n ^ 0. 

k=l,-- ,N, 
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These equations for Qn are inductively integrable. The first (N+1) of solutions are 
given by Qo{t) = Qo{0) and 

ge,(t) = QeM - 1 [ Ad,-.ci^,-.c2(sf) ® )go(o)rfs, k = i,--- ,n. 

Jo 

In order to write down solutions in a compact form, we introduce N linear operators 
Tfc : C°°(R; C^'x^") _ C^IR- C2'x2'"), k = l,--- ,N, defined by 

{%F){t) := -2 r Ad,-.c,^,-.c,(sf ) ® F e C°°(M; C2'>^2'")^ (g^g^ 

Jo 

Then, in terms of these operators, Qe^it), k = 1, ■ ■ ■ , N, take the form 

Q,,{t) = ge,(0) + {TkQomit), k = l,---,N, 

where Qo{0) in the right-hand side is viewed as a constant function. 

We proceed to, say, Qei+e2(t)- The equation for gei+e2(^) is expressed as 

^61+62 — Adg-t4i^g-te2 (S^ ^ ® ^)Qei ~ ^ Adg-t5i^g-te2 (S^ ® ^)ge2) 

and integrated to yield 

= ge,+e2(0) + {T^Q^Mm + {TiQ^MW) + (^2 o Ti Qo(0))(t) + (Ti o Qo(0))(t). 
In a similar manner, the Qn{t) is inductively integrated and expressed as 

where \m\ ■= + ■ ■ ■ + m^, and where Tk : C^'^^"') C°°(M; C^'^^™) are 

defined to be 

Tk := %,o...o r,|^| , = ■ ■ ■ , fc,^,) G {1, ■ ■ ■ , iV}IH. 
Hence, we have obtained the following. 

Proposition 3.1. The solution C{t) to (13. 7p is put in the form of power series in Sk, 
k = lr-- ,N, 

at) = (e*«^ ®e*«^) J] J] J]5"(r;,Q„_^(0))(t) 

n m K 

= (e*«i®e*«^) E e-'iTKCmit), (3.10) 

#{iy|fc^,=fc}=mfc, fe=l,--- ,N 

where we have used the fact that C(0) = J2n^^^n{0) = J^n^^Qni^) ■ 

If the control is piecewise constant in t, the continuation along time of solutions of 
the form f l3.10p with respective constant controls will yield a solution to (13. 7p . 
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3.2 Controllability of the NMR systems 

From Eq. (13.61) . we observe that the control vector field are vertical, so that the measure 
defined in fl2.3l) does not change along this vector field. Hence, the control would not 
make the NMR system entangled without the coupling with the drift vector field. With 
this observation in mind, we gain insight into fl3.10p to understand how the coupling 
occurs. We first look into the operators % defined in (13.91) . Writing out the integrand 
of (13.91) in the form of a power series in s, we obtain 

+ ^ ([6, [ei, )]]F(.)(sf )^ + 2[6, )]F(.) [6, 

+sfV(.)[6,[6,4'¥) + ---}^^- 

By applying Proposition \'2.2\ we can put the above equation in the form 

+ -^{[X^im^ [X^i^i, -X.^{fc)^^(fc)]](F(s)) + 2[X^,^7, -X.^(k)^^(k)]]{F{s)) 
+ [Xj^^,, [Xj^^,, -X^w^^«]](F(s)) )+■■■} ds. (3.11) 

Since the expansion (13.111) means that the drift and control vector fields are coupled 
indeed, and since the operator Tk is the composition of these 7^., Eq. (I3.10p shows 
that the drift and control vector fields work together to give rise to the solution. If 
the coupling between the drift and control vector fields generates vector fields with 
non-vanishing horizontal components, the NMR system will get more entangled by 
controls. 

Incidentally, if the NMR system is controllable, then any initial state can be trans- 
ferred into, say, a maximally entangled state by some controls. This means that the 
control vector field should be coupled with the drift vector field in the case of control- 
lable NMR systems. 

According to [2], the system (13. 3p is said to be pure state controllable if, for any 
initial and final states, there are controls ^q, a = 1, ■ ■ ■ , n, and a finite time T > such 
that the initial state is transferred into the final one on the finite time T. Further, 
the system (13. 3p is pure state controllable if and only if the Lie algebra Cfj generated 
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by {—iHd} U {^crj"'"''}"~]^2" is isomorphic to su(2"'). For the time being, we are 
interested in the couphng between the drift and control vector fields, and postpone the 
controllability problem to Appendix B, in which we will give a necessary and sufficient 
condition for controllability. 

However, it is of great help to apply the controllability theorem to the NMR system. 
To state the theorem, we need to introduce the notion of spin graph. Here, the spin 
graph associated with a given drift Hamiltonian Hd is defined to be the un-oriented 
graph that has n nodes and edges joining the nodes labeled as a and (3 such that 
Jafs 7^ 0, where nodes and edges represent spin-^ particles and interactions among 
spin- 1 particles such that J^/? 7^ 0, respectively. In terms of the spin graph, the 
controllability theorem is stated as follows (see Theorem B.8): The NMR system is 
controllable, if and only if the associated spin graph is connected. In view of this, we 
may set some of Jap in (13.41) to be zero. We here take 

Jipy^O, £ + l<f3<£ + m, 

Ja,e+i7^0, l<a<i, (3.12) 
Jai3 = 0, otherwise. 

Then, the associated spin graph is connected, so that the NMR system with these 
coupling constants is controllable. The drift vector field in (13. 4p is now expressed as 



(l,m)iT 
3 J • 



e+i<i3<e+m 2<a<e 
In this case, one has N = 2, and 

yW _ r^{i.^)i y(^) _ \^ t(i) r^(/3-^).™)iT . t{i) _ j 

e+i<f3<e+m 

2<a<e 

i<a<e e+i<p<e+m 
where ^^,^(3 e span^{iai,ia2}. 

4 Controls in two-qubit NMR systems 

In this section, we deal with the two-qubit case. Our system is expressed as 
dC 
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where J > is the couphng constant between the two qubits, and where ^1,^2 are 
constant controls taking values in spanj^jicri, i(y2}- In comparison to (13.71) . this equation 
is easier to solve, and we can investigate solutions in detail to understand how controls 
make the NMR system more entangled. 

4.1 Two-qubit NMR systems without control 

To understand why one needs controls to make the NMR systems entangled, we look 
into the uncontrolled system with ^1 = ^2 = in (14. ip . If ^1 = ^2 = 0, Eq. (14. II) is 
easily integrated to give the solution 

- (^_,^3^.3(C(U)) - ^ giji/4c^^(o) e-^*/4cn(0) J ' ^^'^^ 
Then the measure F{C(t)) = det (C(t)*C(t)) is evaluated as 

F{C{t)) = |coo(0)cii(0)p + |coi(0)cio(0)p 

- (e-^-'*coo(0)coi(0) cio(0)cn(0)) + e^''*^;MO)coi(0)cio(0)^;i(0))) , 

which gives rise to the inequality 

FiCit)) < (|coo(0)cn(0)| + |coi(0)coi(0)|)2 =: 0(0(0)). (4.3) 

Eq. (14. 3 p implies that the quantity D(C(0)) determines whether the state driven by 
(14. ip with the initial state C(0) reaches a maximally entangled state without any 
controls or not. 

Let us be reminded of the fact that the measure F ranges over < F(C) < 1/4 for 
the two-qubit, where the minimum or the maximum occurs, according to whether C is 
separable or maximally entangled. If D(C(0)) = 1/4, the state C(t) evolves to reach 
a maximally entangled state without control even if C(0) is separable. For instance, 
starting with the separable state 

c(o)=i(; (4.4) 

the solution curve C(t) passes through the maximally entangled state at t = vr/J, 
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We note in addition that F{C{t)) = (1 — cos Jt)/8 in this case. In contrast with this, 
if D{C{0)) = 0, then C(0) is separable, and the state C{t) is always separable because 
of < F{C{t)) < D{C{0)) = 0. An example of such an initial state is 

C{0) = diag (e*^ 0) . (4.5) 

This initial state never gets entangled without controls. From these facts, we observe 
that the separable states i(|0) + |1)) ® (|0) + |1)) and e*^|0) ® |0) are of different nature 
with respect to the drift Hamiltonian. 

4.2 Behaviours of solutions in the two-qubit NMR system 

We turn to the controlled system (14. ip with and .^2 non-vanishing constants. Since 
the associated spin graph is connected, the two-qubit NMR system is controllable, so 
that any separable state can get entangled. This implies that the coupling between the 
drift and control vector fields must occur. 

We solve the Schrodinger equation (14. ip as a power series of the coupling constant 
J > 0. As in the previous section, we assume that the J is small enough, and that the 
solution C(t) of (14. ip can be expanded into a power series in J as C(t) = J2'^=o J^Cnif). 
Substituting this into (14. ip brings about the series of differential equations for C^, 

c'o = iiCo + Cq^2 ■> 

(4.6) 

Cn = -icTsCn-icrJ + ^iC„ + C„^^, n>l. 

These equations can be integrated inductively, like (13. 8p . In fact, from the above 
equations, differential equations for e~*^^C„(t)e~*^2^ are easily obtained and integrated 
to give 

00 

n=0 

= ^^{-iJY / (is„ Ade-s„«i((T3) / Adg-s„_ici (0-3) • ■ ■ / c?si Adg-n^i (as) 

n=0 "^0 Jo 

X C(0) Ad^.^,T(aJ)---Ad^.„_^,.(crJ)Ad^.„,.(4), (4.7) 

where we have used the assumption that C(0) = Yl'^=o J^CniGi). The following propo- 
sition is easy to prove. 

Proposition 4.1. The power series (14. 7p uniformly converges to a function P{t), with 
which the solution of the NMR system (14. ip is given by C{t) = e^^^P(t)e^^^ . 
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So far we have found the solution as the power series (14.71) in the couphng constant 
J. We now analyze the solution in order to see in detail that the coupling between 
the drift and control vectors actually occurs to make the NMR system more entangled. 
We denote the n-th order term in J of (14.71) by Pn(t). Then, the solution is expressed 
as C(t) = J2'^=o'J"^^^^^riit)e^^^ . This explains how the controls and ^2 contribute 
to entanglement promotion of the system. The zeroth-order term does not make the 
system more entangled because of F(e*«iPo(0)e*«2 ) = F{C{0)). Contrarily, the higher- 
order terms in J are expected to make the NMR system entangled in general. We look 
into (14. 7p in detail. The first-order term Pi(t) of (14.71) is put in the form (7iC(0))(t). 
Then, as is seen from (13. lip . Pi(t) is expressed as 



'-lo-3(X>o-3j 



C(0) 



+ ■■■ . 

C(0) 

(4. 



This shows that the drift and control vector fields are coupled indeed. We note that 
[Xg^®/, -Xj^3C3^3]] = [X^i®/, [X/^gj, -Xi^3g3^3]] by the Jacobi identity and that 

By a similar computation, the second-order term P2{t) of (14.70 is written out as 



P2{t) = {T,oT,C{mt) 

^ {-za,{-za,mal)al) - ^3(1^1, cr,]mal + ^3^(0) [6, a,] 



We formally denote the first and second terms of the right-hand side of the above equa- 
tion by ||-(— Xio-3,^o-3)^|c{o) and ||-(— Xjo-3^o-3)([X^j,g,/, —Xi„^(^^.^] + [Xi^^,^, —Xi„^(^^.^])\c(Q), 
respectively. The first term ^{—Xifj^^^^)'^ comes from the Taylor expansion of the flow 
^-ia30<73 with respect to t. 

It is inductively shown that the n-th order term Pn{t) is expanded into a power 
series of t with the lowest-order term ^[—Xi^^^^^Y , which comes from the Taylor 
expansion ^^Ua^f^a^iC {Q)) = J2n=o !iT(-^io3®<x3)"|c(o)- Hence, we have 

Proposition 4.2. If t is sufficiently small, the solution C{t) of (14. ip allows of the 
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approximation of the form 

e-'^^C{t)e-'^^ ~^%^^^{Cm 
Jt^ 



Jt^ ( 



C(0) 



+ [^/(8)6' [-^68>-f' ~-^ifT3®a3]] + [-^/<X)6' [-^/8)C2' ~^ifT3(g)(T3]] ) H ? (4-9) 

/ C(0) 

where ^'^-iazm-.X^^^^ solution of the uncontrolled system dC/dt = —JXi^.^^^^.^ 

with the initial state C(0). 

4.3 Bases of vertical and horizontal subspaces 

So far we have shown that the drift and control vector fields are coupled. We now 
wish to show that the coupling generates vector fields with non-vanishing horizontal 
components. To this end, we have to look into vertical and horizontal vector fields in 
detail. 

Since the subspaces Vc and He depend on the stratum to which C belongs, we 
make a brief review of a stratification of the state space M for the two-qubit system 
after [12]. According to singular values Ai(C) > \2{C) of C G M, the state space M of 
the two-qubit system is stratified into three G- invariant strata as M = Mq U Mi U M2, 
where 

Mo = {C G M I Ai(C) = 1, \2{C) = 0}, 

Ml = G M| Ai(C) > A2(C) > 0}, (4.10) 

M2 = {C eM\ Ai(C) = A2(C) = I/V2}. 
The strata Mq and M2 are the sets of separable and maximally entangled states, re- 
spectively. 

We now deals with bases of the subspaces Vc and He- We take up a typical matrix 
A = diag (Ai, A2) with Ai > A2, to which all the matrices with the same singular values 
are translated by the U(2) ® U(2) action. Let us start with the vertical subspace Va_. 
Proposition 12. II with fl4. lOp implies that 

6, if A G Ml, 
diml^A = <( 5, if A G Mq, (4.11) 
4, if A G M2. 
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Tangent vectors Xu^j{A), Xi^^^j{A) and Xn^^^^A) for j = 1,2,3, defined in 02.101) . 
span the tangent space Va- However, some of these vectors coincide with one another, 

(1) Xi^.^^i{A) = Xij(^^.^{A) at any diagonal matrix A, 

(2) Xi,355/(A) = Xu^aM) = Xij^j{A) at A = diag(l,0) G Mq, 

(3) X,,,55/(A) = Xij^aA^) and X,^^^i{A) = -Xu^^^{A), at A = I/V2 G M2. 
Thus, respective bases of V\ are obtained as follows: 

Proposition 4.3. Let A = diag (Ai, A2). 

(i) If A lies in the principal stratum Mi, then V\ has the basis 

X/(g)/(A), Xjo-3(g)/(A) = Xj/0o-3(A), Xjo-^^/(A), Xj/^o.^(A), Xi^-^^j^A), Xii(g,„^{A). 

(a) If A is a separable state, i.e., A = diag (1,0) G Mq, then V\ has the basis 

Xi/0/(A) = Xjo-3(^/(A) = Xij^„^{A), Xjo-i(g)/(A), Xj/,^o-i(A), Xjo-2(g)/(A), Xn^^^^A). 

(Hi) If A is maximally entangled, i.e., A = I/\/2 G M2, then Va has the basis 

Xii^i{A), Xjo-30/(A) = Xij(g,^^{A), Xjo-j0/(A) = Xj7-0o-j(A), Xi^^siii^) = ^-^j/®o-2(A-)- 

S'mce the G -action (12. 4p isometric, the singular value decomposition (12. 5p . viewed 
as a map, A ^ C = {g ® h) ■ A, gives rise to a basis of Vc from that of Va by the 
differential map {g ® 

The next task is to find a basis of the horizontal subspace He- From (14. lip , the 
dimension of He proves to be 

fl, ifCGMi, 
dimifc=<2, if Cg Mo, (4.12) 
[3, ifCGMa. 

The following proposition is easily verified by a straightforward computation. 
Proposition 4.4. Let A = diag (Ai, A2). 

(i) If A lies in the principal stratum Mi, then the horizontal subspace H\ is a one- 
dimensional vector space spanned by 

Xt7i®cr2(A-) = Xjo-200--^ (A). 

(a) If A is separable, i.e., A G Mq, or A = diag (1,0), then H\ has the basis 

Xiai(^ai{A) , Xj(jiig)cr2 (A) = Xj(j2ig)o-i (A). 
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(in) If A is maximally entangled, i.e., A G M2, or A = I/V2, then H\ has the basis 

(A), iXi^,^i{A). 

A basis of Vc are formed from that of V\ by the differential map {g ®h)^. 

From these propositions, we observe that when A = I / \/2 E M2 = U(2), the 
vertical subspace V\ is identified with the space u(2) of 2 x 2 skew Hermitian matrices, 
and that the horizontal subspace H\ with A = I /V^ is identified with the space isu{2) 
of traceless Hermitian matrices. Further, the basis vectors of H\ given in Proposition 
(iii) are alternatively expressed as 

1.1 1 1 



— ^/) — iXij(^„2{--j=I) — 2Xi„-^^^„.^{--j=I) — 2Xjo-3(g)cri(-^-' 



Remark. The vector fields Xi^^mi ^^'^'^ -^io-i^o-i (11"^ generated, by taking the Lie brack- 
ets among Xifj^^fj^, Xifj-^^j and X^^^^, 

respectively. The first and second equations of the above imply that, if doubly coupled, 
the drift and control vector fields generates a horizontal vector at C = A E Mi and 
another horizontal vector at C = A E Mq, respectively. 

4.4 The concurrence and control 

We now investigate the right-hand side of (14.91) to examine the effect of the control on 
entanglement promotion. We suppose that the initial state C(0) is a diagonal matrix 
A = diag(Ai, A2). 

The terms in the (9( Jt^)-term of the right-hand side of (14.91) are expanded as 

[Xj7-0^2 5 ~-^iu3®a3] = ~X2Xiu^^u^ + ?/2Xjo-3®(Ti , 

where = iXaCi + iyaC2, a = 1, 2, with Xa,ya ^ I^- These vector fields are all vertical 
at A if A 7^ I/V2. Indeed, the tangent vectors Xi^^^xia-A^), Xi„^(^„.j{A), Xi^.^(g,„^{A), 
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and Xjo-3^a-i(A) are written as linear combinations of the vertical vectors Xjo-i(g)/(A), 
Xi^2^/(A), Xij(^^^{A) and Xijf^^^W. For instance, one has 

+ A1A2 
-^i<T2®<T3 (A) = ^^^^^2— — ^Xj^20/(A) + -^Y—j^Xii(g,„^{A), (4.16a) 

A^ + A^ A1A2 

^iai^aai^) = 2(_X2 _ ^2^^«'^l®'f(^) ~ ^2 _ ^2 ^^^^f^l (^) ' (4.16b) 

if A ^ J/v^. If A = I/V2, they are found to be horizontal on account of (14.131) . 
The (9( J^t^)-term is also vertical at A, as is shown by 

1 1 

The (9(Jt^)-term in (14. 9 p is significant for entanglement promotion. The vector 
fields [X^i07, [X^ 
are expressed as 

~^«o-3(g)(T3]] = + ?/^)Xjo-3(gi(T3, 

[-^/(g)^' [^-f<X>6' ~^jo-3®o-3]] = (^2 + y2)^«o-3'»o'3, 



1(T1(X'0'1 ) 



respectively. Since C(0) = A is diagonal, one has Xi„^(g,f^.^{A) = Xii^j{A) G Va, so 
that the tangent vectors [X^^(^i,[X^^(^i, -Xi^.^(^^.^]]\^ and [Xi^^^,[Xi^^^, -Xi^^^^^]]\^ 
are vertical. In contrast with this, from Proposition l4.4l and Remark after it, it turns out 
that the tangent vector [X^^^j, [Xj^^^, — -^io-3(g.(T3]] |^ contains the horizontal components 
-^jo-2®o-i(A), Xjo-i(g)o-2(A)- In particular, if and ^2 are taken as ^1 = iai and ^2 = 
respectively, the tangent vector in question is horizontal, 

[XicTi^I, [XiJ^fj^, — Xjo-3(g)(73]] = Xia-2^ai{A) E H\. 

This shows that the drift vector field Xi„^(^(j^ and the control vector fields X^-^^j, Xj^^^ 
are coupled to generate a horizontal vector at A. We note further that the horizontal 
vector emerges at the third-order term in t, if A 7^ I/V^- 

We wish to evaluate the entanglement to confirm that the two-qubit system gets 
more entangled. As our entanglement measure and the concurrence are equivalent 



for the two-qubit system on account of -\/ F{C{t)) = \ det C{t)\, we here use the con- 
currence as a measure. From Eq. (14. 9p . the concurrence is approximately evaluated 

as 

^MCW) - ^jF{^il^^^^{A)) + ^^J^^ {dF)A{[X^,^j, -X..3^<x3]]) , (4.17) 
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where we have used the following facts; (i) F{e ^^^C(t)e ) = F(C(t)), (ii) [X^^^j, -X, 
and — Xjo-300-3] are vertical at A if A is not the maximally entangled state 

I / \/2^ (iii) — Xjo-3®(73 i^^iiX)/) ""^icrsixxTs]] and —Xi^^^f^^\Xji2)^^^—Xi^^^„^ are vertical at 
A, (iv) [X^^(^i,[X^^^i,-Xi„.^^„^]] and [Xi^^^,[Xi^^^, -Xi„^^„^]] are vertical at A, and 
(v) [Xgi®/, -Xi^355^3]] = [X^i^/, [X/cg^^, -Xj^g^^g]]. Since the derivative of the 

measure F is given by 

dF = 8Re (det C* tr (cxsCasrfC^) ) , (4. 18) 

one has 

{dF)A([X(i(^i, — Xj^ggi^g]]^ = — il/2^ 1/1 _ _\2^_ 

Hence, substituting this into (14.171) . we have the following estimate of the concurrence, 



^^FiCW) - ^FiA^^ZjM) + + yix,){\l - \l) (4.19) 

for a sufficiently small t > 0. It is to be noted that the effect of controls on the 
concurrence emerges actually at the third-order term in t. We have to note that for 
A = diag(l,0), one has F{ipi\^,^^^.^{K)) = 0, so that ^F{C{t)) ~ ^{-xm + x^yi). 
This means that the state A = diag (1,0) gets entangled slowly by the controls. 

The growth rate of the concurrence at the third order in t depends on the quantity 
—Xiy2 + yiX2, where Xa.yp are control parameters given by = ix^cri + iyaO'2- It 
would be reasonable to restrict the magnitude of controls to Yla=ii^a + = 1- If 
this is the case, the maximal growth ratio is realized when xi = —y2, X2 = yi- This 
means that the most efficient control for entanglement is given by ^1 = ixiai + 1X2^2-, 
^2 = ix2<yi — ixia2 with x\ + x\ = 1/2. However, if we choose the control given by 
^1 = ixiai — ix2<J2, C,2 = ix2<Ji + ixia2 with x\ + X2 = 1/2, the two-qubit gets into a 
less entangled state. If ^1 = ixiai + ixia2, ^2 = ix2(Ti + 1x2(72 with xf + x^ = 1/2, no 
change will occur in the entanglement of the two-qubit. 

We turn to the case that the initial state C(0) is not diagonal. Let us be reminded of 
the fact that when the C(0) is diagonal, the (9(Jt^)-term in (14. 9 p is vertical. However, 
if C(0) is not diagonal, the (9(Jt^)-term is not vertical, so that Eqs. (I4.15P and (I4.18P 
are put together to yield the contribution to entanglement promotion by 



2,/F(cW) 



{dF)c{0) (^X^^^I, — Xj^^g^^g] + [X70^2, — Xj<ygg,(^g]^ 



- ^Re( e-'%coi + cio)(ircoo + ^2011) ) , (4.20) 
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where = Xa + iya for a = 1,2, and where e*^ is defined by detC = e*^| detC| = 
e^^ \JF{C). From fl4.20p . the (9(Jt^)-term identically vanishes for arbitrary and .^2 if 
Coo = Cii = or coi + Cio = 0, but it takes a non-zero value in general, if Cqi + Cio 7^ 
and (coo,cii) 7^ (0,0). 

5 Concluding remarks 

We have discussed how the control and drift vector fields are coupled to make the NMR 
system more entangled on the two-qubit model, by examining the solution in detail in 
terms of horizontal and vertical vector fields and by evaluating the concurrence. As 
was pointed out in Sec. 4.1, the initial state C(0) = diag (e'^,0) or e*^|0) © |0) never 
gets entangled without controls, but we have shown in Sec. 4.4 that it can be entangled 
if controls are taken. However, the effect of the control on entanglement is so slow that 
it emerges at the third-order term in t. In fact, as is shown in Sec. 4.4, only vertical 
vector fields are generated in the second-order term in t by the coupling between the 
drift and control vector fields, and further vector fields with non- vanishing horizontal 
components emerge at the third-order term in t. 

It would be of help to state the relation of our theory to the Cartan decomposition 
of the Lie algebra. As for the two-qubit systems, the associated Lie algebra is su(4), 
of which the Cartan decomposition is given by 

5u(4) = ! © p, 

where 

t = span{zJ © aj,iak © I}j,k=i,2,3, P = span{2iaj © ak}j,k=i,2,3- 

The subalgebra 6 generates vertical tangent vectors, and some of elements in p span the 
horizontal subspace at A = diag (Ai, A2). For example, if Ai 7^ A2, one has a horizontal 
vector A'jo-i(X)cr2(A) (see Props. 4.3 and 4.4). However, the vectors associated with p are 
not always horizontal. As is seen from (14.161) . Xi^^t^^j.^^A) and Xi^^^i^^^^A) are vertical 
if A ^ I/y/2. In [27], as to non-local operation, they state that t can be viewed as the 
local part in 3u(4) and p as the non-local part. In this sense, the horizontal vectors are 
of non-local nature, and the horizontal subspace makes the non-locality quite sharp. 
This is because the entanglement gets promoted most efficiently in the direction of 
horizontal vectors. 
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We comment on the motivation behind this study, which is related to the decompo- 
sition method for the analysis of quantum control problems on compact Lie groups. The 
decomposition method is of much use to design a control in order to construct a desired 
unitary operator [15-17], and then has been developed [7,8]. In these situations, the 
control (or local operation) works instantaneously, and the drift (or non-local) Hamil- 
tonian is assumed to be negligible while the control works. Under this assumption, the 
time evolution operator Ut = e~^^^ can be decomposed into the product of the local 
and non-local operators, 

Ut = Lo e-'^'^" Li • • • e"^"^^'^ Ln (5.1) 

where Lq, ■ ■ ■ .L^ G SU(2^) ® SU(2™) are generated by the control Hamiltonian, and 
where ri, ■ ■ ■ ,r7v > are time intervals for which the control Hamiltonian vanishes. 
In this view, the bipartite entanglement is not concerned by the coupling between the 
drift and control Hamiltonians in the time evolution. Our motivation is to consider the 
coupling in the evolution without the above assumption, and to investigate the relation 
of the coupling to the bipartite entanglement. The work [24] by Romano shares the 
motivation with us. 

In [24], the decomposition of the unitary group SU(4) is used to study entanglement 
magnification. The time evolution Ut is decomposed into Ut = LtAtKt, where Lt,Kt G 
SU(2) ® SU(2) and At E H with H the Cartan subgroup of SU(4). The author says 
that Lt is irrelevant, and claims that the non-local operator At makes a contribution 
to entanglement magnification together with a help of Kt. The decomposition method 
is sharp to study the entanglement in the evolution, because the evolution operator 
is completely separated into the "entangling part" At and the two local operators Lt 
and Kt. But, the process of entanglement promotion does not explicitly appear in the 
decomposition, since not only At but also Lt, Kt take in the coupling between the drift 
and control Hamiltonians. In fact, the At, Lt and Kt are all unknown until the time 
evolution operator Ut is integrated in an explicit form. 

On the other hand, the method of this article is represented in a familiar form with 
bipartite entanglement, and explains how the coupling between the drift and control 
Hamiltonians occurs in the time evolution. See fl3.10p . (13.111) . (14. 7p and (14. 9p . In 
particular, the right hand-side of (14. 9 p is not discussed in the context of (15. ip . In 
a comparison to [24], our method allows one not to solve the time evolution of the 
system, although it is not sharp in the sense that the terms in (14.90 contain vertical 
(or local) components. 
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The formulas, fl3.10l) together with fl3.1ip . and fl4.9l) also show how the effect of 
the controls spreads among qubits in the dynamical system. In contrast with this, the 
proof of the controllability of the NMR system would show the spread of the controls 
among qubits in a "static" situation. For the sake of self-containedness, we give the 
proof of the controllability of our NMR system (see Theorem lB.81) in a method different 
from that in [1]. As is pointed out in App. B, the NMR system is controllable if and 
only if the associated spin graph is connected. 
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A Lie brackets in the two-qubit NMR system 

A straightforward calculation along with Corollary 12.31 provides 
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and the other Lie brackets prove to vanish. Further, double Lie brackets are given as 
follows: 
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B Controllability of the NMR system 

As stated in Sec.3, we study the controllability of our NMR system. The drift Hamil- 
tonian we treat is given by 



Hd = Ho = ^2 Jal3 O's" 0": 
l<a<(}<n 



3 ; 



(B.l) 



where, for simplicity, we have denoted the drift Hamiltonian by Hq and abbreviated 
crj"'"'^ to where aj"'"''' were defined in (13. 2p . For our convenience, we express the 
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control Hamiltonian as follows: 

n 2n 



fi=l fj,=n+l 



where 



(p) _ 1 

H ■= ' ^~ (B 2) 



Our control system with the Hamiltonian H = + He is (pure state) controllable if 
and only if the Lie algebra generated by {Hq, iH^, /i = 1, ■ ■ ■ , 2n, is equal to su(2"'). 

For comparison's sake, we here quote the drift and control Hamiltonians studied in 
[1], which are given by 

a</3 

and 

n 

^i = E""^f' ^■ = l'2,3, (B.4) 

a=l 

respectively, where Map, Nap, Pap £ are coupling constants and G M are the 
gyromagnetic ratio of the a-th spin-| particle (or qubit). In our case of ( ]B.2p . it is 
assumed that each of spin-| particles can be stimulated by the only two components of 
the magnetic field. In contrast with this, in the case of ( ]B.4p . all the spin-| particles are 
assumed to be stimulated simultaneously by all the three components of the magnetic 
field with possibly different gyromagnetic ratios. While our drift Hamiltonian ( ]B.1|) 
contains only (T3 factors, but theirs ( ]B.3p has all the (Tj factors. In [1], they prove that 
the NMR system with the Hamiltonian H'q + X]j=i ^i-^j controllable if and only if 
the associated spin graph is connected. 

Our objective in what follows is to show that our NMR system is controllable if 
and only if the associated spin graph is connected. The main point is to show that the 
Lie algebra generated by the operators —iH^, /i = 0, 1, ■ ■ ■ , 2n, is equal to su(2"). In 
this respect, what we have to do for proof is the same as that in [1]. However, as our 
Hamiltonian is different from that in [1], the method for proof should be different from 
that in [1]. 

B.l A decomposition of 5u(2") 

Before calculating commutators among {— 'j-f^^}^=o,i,- - ,2n, we decompose su(2") into 
the sum of subspaces. We denote by B the Lie algebra su(2"), 

B = sp&n^{iaj, ® a,, ® • ■ ■ ® cr,„ | ji, J2, ■ ■ ■ , Jn e {0, 1, 2, 3}} \ {i/®"}, (B.5) 
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and break it up into subspaces B^^\ k = 0,1, ■ ■ ■ , n — 1, all the basis elements of which 
have as many as k /-factors. Since dimS'-'^-* = (^)3"'~'^, and since dimsu(4) = 4" — 1 = 
YlkZo (fc)'^""'^' vector space B is broken up into 

n-l 

i3 = 0i3('^\ (B.6) 

k=0 

For k = 0, B^^'' is expressed as 

^(0) = span^iicxj, (g) cxj^ (g) ■ ■ ■ (g) aj„ \ ji, j2, ■ ■ ■ , jn G {1, 2, 3}}. 
Among B^'^^ with A; = 1, ■ ■ ■ , n, we write down B^''^ only for k = n — 1; 

= spe.n^{zaf\zaf, zaf \ j = 1, 2, 3}. 

Further, we set 

B3 := span^{iaj, ® (x,, ® ■ ■ ■ ® | Ji, J2, ■ ■ ■ , Jn G {0, 3}} \ {2/®"}. (B.7) 

The following lemma shows that B^ holds a key position in calculating commutators. 

Lemma B.l. Commutators among {— ii/^};x=i,... ,2nUi33 generate all the basis operators 
ofB. 

Proof. Note that we have B3 = T.k=oi^s ^ ^^''^) ^0"^ (H) and flRTl) . For A; = 0, we 
take an operator icrf" G ^3 fl S*-^-' . The commutators between iaf" and —iH^ are given 
by 

. _, ^ 1 _ ® -2 ® -s^^"-'^^ = 1, ■ • ■ , 

^ ^' ^^"\-za3«('^— ^)®ai®a3®('"-'^\ ^ = n + l,---,2n. 
Taking successive commutators [[^cr®", —iH^], —iH^], /x, z/ = 1, ■ ■ ■ , 2n, and so on, we 
can obtain the operators of the form, up to ± sign, 

iaj^ (g) aj^ (g) ■ ■ ■ (g aj„, ji, j2, • ■ ■ , jn G {1, 2, 3}. 

If we start with il (g a®^" G -B3 n B^^^ and follow the same procedure as above, we 
can obtain 

il (g(Tj^0---(g aj^, j2, ■ ■ ■ , in G {1, 2, 3}. 

Through the process of taking the commutator of an operator in B3 fl B^'''^ with —iH^, 
the resultant tensor product operator may have the /-factors fixed and cr3-factors in 
the site yU or fi — n changed to ai or o"2, according to whether fi G {1, ■ ■ ■ ,n}, or 
/i G {n + 1, - ■ ■ , 2n}. Taking successive commutators, we can obtain the tensor product 
operators with the {n — k) cr3-factors replaced by ai or a2 and the k /-factors fixed, 
so that we have all the basis elements of B^''^ in the form of commutators among 
{-i^^}^=i,...,2n U {B3 n B^^^), = 0, 1, ■ ■ ■ ,n - 1. This ends the proof. □ 
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B.2 Controllability with the complete spin graph 

We now assume that J^^ 7^ for all 1 < a < /? < n, or that the associated spin graph is 
complete, to show that the commutators among {— i-ff^}^=o,i- - ,2n span su(2"). On ac- 
count of Lemma FB. 11 we have only to show that commutators among {— 'j-ff^}^=o,i,- -,2n 
can span B^. We first show that by taking commutators among {—iHij,}^=Qx---,2n, the 
term —iaj^^'^o'll^^ of the drift Hamiltonian —iHo can be singled out, if Jaj3 7^ 0. For 
—iHo = —i X]a</3 Jaf3(y^i^ (y^^ and — i-ff^ = \ n = 1, - ■ ■ ,n — 1, the commutators 
between them are given by 

At— 1 n 

0=1 /3=Ai+l 

We calculate further commutators to obtain, for 1 < ^ < u < n, 

[[-iHo, -iH^i], -iHv] = -iJ ^.v'y'i''' <^2^ ■ 
Moreover, a calculation with the right-hand side of the above provides, for 1 < /i < 

If J^v 7^ 0, we obtain —ia^^^^aj^^ as a commutator among {iHo, ■ ■ ■ ,iH2n}- Thus, we 
have proved the following 

Lemma B.2. If Jap 7^ for a and (3 with 1 < a < (3 < n, the operator —ia^^a^'^ can 
he realized as a commutator among {—iH^}fj,=o,i,---,2n- 

We now show that the B3 can be generated from {—iH^}^=o,i,--- ,2n, by using Lemma 
IB .21 under the assumption that Jap 7^ for all 1 < a < /3 < n. Let ai, i = 1, ■ ■ ■ ,m, 
be positive integers such that 

1 < tti < ^2 < ■ • ■ < ttm-i < «m < n with m > 2. (B-8) 
For a given sequence of ctj. Lemma FB.21 provides us with the operators 

in the form of commutators among {—iH ^} ^=o,i,- - ,2n- We can form commutators among 
these operators and {—iH^}ij.=o,i,- - ,2n to obtain 

P ■— r •TT ] _ ■iai) (02) 

T? ■— W ^^i'^k)Aak+i) „-/V 1 ■TT 1 _ „• (Q!fe)^{"fc+i) 

^^•--[[-^0-3 ^3 ,-tHaJ,-lHn+ak+i\ - -iCr2 ^1 

rp r ■ _(Q!m-l) _(am) • rV ] _ -(arn-i) (am) 
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where = 2, ■ ■ ■ , m — 2, if m > 3. From Fi, F2, ■ ■ ■ , -Fm-i, we can form commutators, 
for m = 3, ■ ■ ■ , n, 

Cm :=[■■■ [[Fi, F2], ■ ■ ■ , Fm-2], Fm-i] = ■ ■ ■ 

which is in B3 n For m = 1, 2, the operators Ci, C2 can be expressed, respec- 

tively, as 

Ci := = l<a<n, (B.9) 

C2 := (B.IO) 

where Ci G ^BaH^B^'^-^) and C2 G SsH^B^'^-^), and these can be expressed as commutators 
among {-zi^^}^=o,i,.-,2n as well. 

Taking all possible sequences of ai subject to ( IB. 81) . we can construct all possible 
tensor product operators with m (Js-factors and {n — m) /-factors as commutators 
among {—iH^}^=Q^i^... 2n- In other words, we can form all the basis elements of ft 
g(n-m) jj-^ ^Yie form of commutators among {— ii//i}/i=o,i,---,2n- Since = J2m=ii^3 ^ 
g(n-m)^^ we obtain the following 

Lemma B.3. // Ja/s 7^ for all a and j3 with 1 < a < (3 < n, the B3 are generated by 
commutators among {—iH^}ij.=o,i,---,2n- 

From Lemmas IB. II and IB. 31 it turns out that if Jap 7^ for all a and P with 

l<a<l3<n, Bis generated by taking commutators among {— «//^}^=o,i,-- ,2n- 
Hence, we obtain 

Proposition B.4. // J^js 7^ for all a and (3 with 1 < a < (3 < n, or if the spin graph 
associated with the drift Hamiltonian is complete, the Lie algebra su(2") is generated 
from {-zi/^}^=o,i,-,2n- 

B.3 Controllability with a connected spin graph 

In the following, we consider the case where some of Jo/3's may vanish, but the spin 
graph is connected. Suppose that Jq/3, Jp^ 7^ 0, and Ja-y = 0. We may assume that 
a < /5 < 7. By Lemma lB.21 we can put the operators —ia^'^a^\ —ia^^'^a^^ in the 
form of commutators among {— 'ji/|j}^=o,i,---,2n- We then take commutators to obtain 

D := -^H,], -^af4^^] = (B.ll) 
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From D, —ia\ Vg , —iHn+p^ and —iH^, we obtain 

D' := A[[D, -zH,i [-iH^^p, -^4^^^^^]] = -laf^ a^^\ (B.12) 

Further calculation provides 

D":= [D',-iH^^,] = ~iat^a^^\ 

Thus we have found that the operators —ia^^'a^^^ can be described as commutators 
taken among —i(j^''a^\ —ia^f''a^\ and —iHfj,, = 1, ■ ■ ■ ,2n. This implies that if 
there are interactions between particles a and (3 and between particles (3 and 7, an 
interaction between particles a and 7 is induced, though Ja-y = at the beginning, by 
taking commutators among {— i-ff^}^=o,i,--- ,2n- We may interpret this fact as follows: In 
the process of making commutators among {— «i^^};i=o,i, --,2n5 two edges 7) 
of the spin graph give rise to a new edge {a, 7) which represents an interaction between 
two particles a and 7. 

So far we have studied the case where there are two edges between the nodes 
a and 7. We may apply the same procedure to the case where there are r edges 
between the nodes a and 7, where r = 2, ■ ■ ■ ,n — l. Suppose we have, say, three edges 
{a, Pi), {^1,^2), {^2,1) with Jai3i, <//3i/32? ^1321 0- The same procedure as above yields 
the edge (a,/52) from {a, Pi) and {Pi,P2), and consequently (a, 7) from {a,P2) and 
{P2,'j)- Then, it turns out that if two particles a and 7 are linked with a sequence of 
two-particle interactions, the operator —icx^^^a^^^ describing an interaction between the 
particles a and 7 is induced by taking commutators among {— ^-ff^}^=o,i,- - ,2n- Since 
the spin graph is connected, there is a sequence of edges between any pair of nodes 
a and 7, so that one can obtain the operator — icTg^Vg^'' expressed as a commutator 
among {— 'i-ff^}^=o,i,- -,2n- Thus we have reached the same conclusion as in Lemma [6.21 
with an assumption weaker than that in Lemma IB.2[ 

Lemma B.5. // the associated spin graph is connected, all the interaction operators 

n, are generated from the drift and control Hamiltonians. 

Thus Prop. IB. 41 is refined as follows: 

Proposition B.6. If the spin graph associated with the drift Hamiltonian is connected, 
the Lie algebra su{2'^) is generated from {— 'i-f^^}^=o,i,---,2n- Thus, the NMR system is 
controllable if the associated spin graph is connected. 
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B.4 A necessary and sufficient condition for controllability 

We now consider the case where the spin graph S is disconnected. We assume that 
the graph S is broken up into two disjoint subgraphs, Si and S2, each of which is 
connected. Suppose that Si and S2 have nodes 1 to r, and r + 1 to n, respectively. 
According to the decomposition, S* = 5*1 U 6*2, of the spin graph, the drift Hamihonian 
is also decomposed into the sum of two terms, 

l<a<l3<r r+l<a<l3<n 

Since 5*1 is connected, in the same procedure as taken in proving Prop. IB. 61 the Lie al- 
gebra Ci generated from —iH'Q and {—iHa, —iHn+a}a=i,---,r proves to be Ci = su(2'') ® 
j^{n-r)^ In the same manner, the operators —iHQ and {—iHr+a,—'>'Hn^r+a}a=i,---,n-r 
generate the Lie algebra C2 = ® su(2"'~^'). Since two operators each of which is 
in the respective Lie algebras £1 and £2 commute, the Lie algebra generated by the 
whole operators {—iHfj.}^=o,i,-,2n splits into the direct sum Ci © £2, which is a subal- 
gebra of su(2"') but not equal to the whole su(2"'). We may generalize this fact to the 
case where the spin graph is broken up into more than two disjoint subgraphs, but the 
generalization is easy to perform. So far we have shown the following 

Proposition B.7. // the spin graph associated with the drift Hamiltonian is discon- 
nected, the NMR system f l3.3l) is not controllable. 

From Props IB. 61 and IB. 71 we obtain the following theorem; 

Theorem B.8. The NMR system (13.31) with the drift and control Hamiltonians (13. ip 
is controllable, if and only if the spin graph associated with the drift Hamiltonian is 
connected. 



